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New physics from ultrahigh energy cosmic rays 
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Cosmic rays from outer space enter the atmosphere with energies of up to 10 GeV. The initial 
particle or a secondary hadron inside the shower may then interact with an air nucleon to produce 
nonstandard particles. In this article we study the production of new physics by high energy cosmic 
rays, focusing on the long-lived gluino of split-SUSY models and a WIMP working as dark matter. 
We first deduce the total flux of hadron events at any depth in the atmosphere, showing that 
secondary hadrons can not be neglected. Then we use these results to find the flux of gluinos 
and WIMPs that reach the ground after being produced inside air showers. We also evaluate the 
probability of producing these exotic particles in a single proton shower of ultrahigh energy. Finally 
we discuss the possible signal in current and projected experiments. While the tiny flux of WIMPs 
does not seem to have any phenomenological consequences, we show that the gluinos could modify 
substantially the profile of a small fraction of extensive air showers. In particular, they could produce 
a distinct signal observable at AUGER in showers of large zenith angle. 

PACS numbers: 96.50.sd, 12.60.-i 



I. INTRODUCTION 

The standard model (SM) has been very succesful 
when confronted with the experimental data from parti- 
cle colliders. Despite that, not much is known about the 
nature of the Higgs sector. Hopefully, the Large Hadron 
Collider (LHC) at CERN will provide the value of the 
Higgs mass and will unveil the mechanism that explains 
that value. A dynamical mechanism would imply new 
physics just above the electroweak scale: SUSY, little 
Higgs, technicolor, extra dimensions, or some other un- 
expected physics. It is also possible, however, that the 
LHC just completes and confirms the SM as it was formu- 
lated, pushing above the TeV scale the limits for this new 
physics [l], Q. This possibility has been seriously con- 
sidered after recent astrophysical and cosmological data 
suggested a non-zero vacuum energy density [3|. The 
preferred value A w (10 -3 eV) 4 does not seem to be ex- 
plained by any dynamics at that scale, and could be just 
telling us that the universe is larger and/or older than 
we expected 1 : a multiverse Q. From the experimen- 



tal point of view, precision data (flavor physics, EDMs) 
[1] suggest that the scale of new physics is well above 
the TeV, whereas neutrino oscillations point to an even 
larger scale, above 10 11 GeV [6|. In contrast, dark mat- 
ter strongly suggests a stable WIMP below (or at) the 
TeV Therefore, the question of what (if any) physics 
beyond the SM we should expect at accessible energies is 
wide open. 

On the other hand, it is known that accelerators are 
not the only place to look for collisions of TeV energy. 
There is a well established spectrum of cosmic rays that 
extends up to E — 10 11 GeV [8|. When a proton from 
outer space enters the atmosphere and hits a nucleon, 
the center of mass energy i/s = \f2m^E in the collision 
may go up to 500 TeV. Depending on the energy of the 
initial particle, after the first interaction there will be 
many secondary hadrons 0, with still enough energy 
to produce TeV physics in a collision with an atmospheric 
nucleon. If created, the massive particle would be inside 
the air shower together with thousands of other standard 
particles. Therefore, to have any chance to detect it the 
new particle should be long lived (its decay products in- 
side the shower would be undetectable). This is precisely 
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A could take different values in different regions. To expect one 
region (the universe we see) with a value 10 120 times smaller 
than its natural value, around Afp ;aJlcfc , there should be over 



10 120 of these regions. Analogously, our universe could have been 
preceded by 10 120 Big Bang/Big Crunches of different vacuum 
energy. 



2 



the case of the gluino of split-SUSY models Q. Such a 
particle could change substantially the profile of the air 
shower where it was produced or it may survive (together 
with muons and neutrinos) the shower and be observed 
far from the initial interaction point, in quasi- horizontal 
events or in neutrino telescopes. 

In this paper we study the production of new physics 
by cosmic rays. To evaluate that we determine the to- 
tal flux of hadrons in the atmosphere: primary plus sec- 
ondary nucleons, pions and kaons. The decay length of 
these particles at these energies is much larger than their 
interaction length in air, so the probability that they pro- 
duce new physics is just 

V X = , (1) 

where ax is the cross section to produce the exotic par- 
ticle and <tt the total cross section of the hadron with 
the air. We deduce as well the spectrum of secondary 
hadrons created by a single proton of ultrahigh energy. 
We then show how to apply these results to evaluate the 
production rate of long-lived gluinos and also of a WIMP 
working as dark matter of the universe, and discuss pos- 
sible signals in air showers and neutrino telescopes. 

Other studies of the production of new physics by cos- 
mic rays refer to primary neutrinos [111 ]. In particular, 
there are several recent analyses of the possibility to ob- 
serve long-lived staus at neutrino telescopes [12| . A neu- 
trino of very high energy may reach a telescope and in- 
teract there with a nucleon. Since the neutrino is weakly 
interacting, the relative effect of new physics is here more 
important than for protons or electrons. If IceCube 13 1 
establishes the neutrino flux at these energies, the search 
for deviations to the standard neutrino-nucleon cross sec- 
tions will put strong constraints on some types of TeV 
physics [141 ]. 



II. TOTAL FLUX OF HADRONS 

The flux of nucleons reaching the atmosphere in the 
range of energy from several GeV to ss 10 6 GeV can be 
approximated by 0] 

d&N - _ „_o •> nucleons 



10 



1!) 



1.8 E- 2 - 1 



(2) 



dE cm 2 s sr GeV ' 

with E given in GeV. Around 90% of these primary nu- 
cleons are protons (free or bound in nuclei) and 10% neu- 
trons. At energies « 10 6 GeV (the knee) the differential 
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FIG. 1: Flux of primary nucleons (solid), total nucleons, and 
secondary pions and kaons. 

spectral index changes from 2.7 to 3, and at 10 10 GeV 
(the ankle) it goes back to 2.7 (see Fig. [T]). 

As this flux of nucleons enters the atmosphere, its in- 
teractions with the air will induce a flux of secondary 
and tertiary hadrons (all of them referred as secondary) 
at different vertical depths. Although approximate an- 
alytic expressions can be obtained in limited regions of 
energy 0, H3] , it is easy as well to perform a numerical 
simulation. 



We use C0RSIKA 



151 ] to obtain the number and energy 



distribution of secondary nucleons, (charged) pions and 
kaons (we neglect the baryons containing an s quark) pro- 
duced inside the air shower started by a nucleon of energy 
E. In Fig. [T] we plot the total (primary plus secondary) 
flux of hadrons with energies from 10 4 to 10 11 GeV. We 
obtain this flux generating 1500 showers per decade of 
energy. We observe that at energies around 10 7 GeV 
(t/s i=s 5 TeV) the secondary nucleons increase in a 50% 
the number of primaries. Note also that the number of 
mesons (all of them secondary) is approximately a 15% 
the number of primary nucleons of the same energy. 

The region of energies beyond 10 8 GeV (the one rel- 
evant at AUGER [l6|) deserves a more detailed anal- 
ysis. We notice an important factor that distinguishes 
two types of showers: the presence or not of a leading 
hadron carrying a significant fraction of energy after each 
hadronic interaction. In Table [I] we give two examples of 
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Energy [GeV] 


N 
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K 
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10 4 -2.1 x 10 4 


7539 


54585 


10092 


9873 


70692 
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2.1 x 10 4 -4.6 x 10 4 


3907 
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4.6 x 10 4 -10 5 


2015 
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3360 
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6771 


1323 
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1634 


2.1 x 10 5 -4.6 x 10 5 


469 


3352 
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4.6 x 10 5 -10 6 


254 
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300 
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10 6 -2.1 x 10 6 
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TABLE I: Spectrum of hadrons in a 10 10 GeV shower with a 
first interacion of elasticity xf ~ 0.3 (left, shower (a) in the 
text) and xf ~ 0.03 (right, shower (bj). 

showers, both of energy E = 10 10 GeV, with a very differ- 
ent spectrum of secondary hadrons. Clearly, shower (b) 
does not include a leading hadron after the first interac- 
tion. A relevant parameter to characterize this feature 
is the elasticity xp of the first hadronic interaction: in 
example (a) we have xp ss 0.3 {i.e., a secondary hadron 
keeps 30% of the initial energy) whereas (b ) corresponds 
to x F ~ 0.03. The analysis of 500 showers of 10 10 GeV 
shows that around 15% of them have xf < 0.06 (i.e., no 
leading hadron after the first interaction with the air). 

III. FLUX OF MASSIVE LONG-LIVED 
PARTICLES 

Let us now estimate the production rate of new physics 
in the collision of these cosmic hadrons with atmospheric 
nucleons. The probability that a hadron h of energy E 
produces the exotic particle(s) X is determined by the 
cross section (E) and the total cross section with the 
air, a% a {E): 

A rr hN 

K(E) « ^g- , (3) 

where we assume A = 14.6 nucleons in a nucleus of air 
and neglect nuclear effects in the interaction to produce 
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C h [mb] 


C\ [mb] 


C 2 h [mb] 
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185.7 


13.3 


0.08 


TV 


100.5 


16.9 


0.00 
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79.7 


13.9 


0.05 



TABLE II: Constants defining the total cross section with 
the air. 
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FIG. 2: Interaction length in air Xq of the gluino hadron. 
The interaction length in water is a 3.2% shorter. 

X. Taking the total fluxes d^h/dE of the three species 
of hadrons in Fig. [T] we obtain that the flux of exotic 
particles (of any energy) is just 

E / dE^-V x (E). (4) 

If we are interested in the energy distribution of X we 
must evaluate 

where Ex is the energy of X and the cross sections in- 
side the integral are evaluated at E. The (default) cross 
sections with the air used by C0RSIKA above 10 4 GeV 
can be approximated by at^ « C# + C\ log(E/GeV) + 
C£log 2 (£:/GeV), with the constants given in TablcHT] 

As we mentioned above, to be of interest the exotic 
particle should be 'penetrating and survive the air shower 
where it was produced. We will then consider two differ- 
ent cases. The first one is the long-lived gluino g of split- 
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FIG. 3: Total cross sections ax N with X = gg and XX f° r 
h = N,tt,K. 
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FIG. 4: Flux of gluino pairs as a function of the gluino mass. 
We separate the contribution of the different hadrons to this 
flux 



SUSY models 0] • In the most evasive scenario the gluino 
fragments always into an electrically neutral gluino-gluon 
hadronic state G (a particular type of R-hadron). The 
experimental bounds on its mass are in this case around 
170 GeV [ljj. Notice that general bounds [18j or bounds 
based on the observed delay in the time of flight versus 
a muon [l9| do not apply to a neutral long-lived hadron. 
The gluino hadron will interact often with the air, its 
interaction length (in Fig. [2|) can be estimated [20| as 
Xq w (16/9) Att at Eq w E^Mjm^ (i.e., when both 
hadrons have the same velocity). However, in each in- 
teraction it will lose a very small fraction of its energy 
AE/E w k/M, where k w 0.14-0.35 GeV and M is the 
gluino mass [ljj, |20j, [21| . Therefore, the gluino hadron is 



very penetrating. In particular, it keeps going after the 
atmosphere has absorbed most of the hadrons in quasi- 
horizontal showers. 

The parton-parton cross sections to produce gluino 
pairs can be easily obtained taking the limit of large 
squark mass in the general expressions given in Ref . 22| , 
and we include them in the Appendix. In Fig. [3] we plot 
the total hN (h = N, it, K) cross sections a-- 1 to produce 
gluino pairs for values of the hadron energy between 10 4 
and 10 11 GeV and M = 200, 300 GeV. We have as- 
sumed isospin symmetry to deduce the kaon PDFs and 
have taken a renormalization scale fj, = 0.2M, as sug- 
gested by a NLO calculation 23 1. 



We can now calculate the flux of gluino pairs pro- 
duced in the atmosphere. In Fig. 2] we plot the flux 
for values of M between 150 and 450 GeV. We express 
this flux in gluino pairs per year, squared kilometer and 
sterad. Multiplying by 27r sterad, for example, we can 
read that the « km 2 IceCube detector would be exposed 
to one downgoing gluino pair per year if M ~ 160 GeV. 
Around a 64% of the gluino flux would be produced by 
the primary nucleons, whereas the rest corresponds to 
secondary nucleons (16%), pions (16%) and kaons (4%). 
In Fig. [5] we give the energy distribution of the gluinos 
for M = 200, 300 GeV. The energy E is in this case 
the total energy of the two gluinos in the lab frame. In 
the next section we comment on the phenomenological 
relevance of these results. 

The second nonstandard particle that we would like to 
consider is a stable WIMP \ constituting the dark mat- 
ter of the universe (see [24| for a recent review on dark 
matter). To define a simple scenario, we assume that \ 
is a Majorana fermion with only weak interactions and a 
relatively light mass, M x « 70 GeV, and that all other 
possible particles (neutral or charged in the same SU (2) i 
multiplet) are at least 30 GeV heavier. The later require- 
ment is necessary to avoid collider bounds, and will imply 
that coannihilations at temperatures below M x are irrel- 
evant. We also need a discrete symmetry similar to the 
R-parity to make x stable. Under these hypotheses, the 
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FIG. 5: Energy distribution of the gluino and WIMP events. 



only relevant parameter to determine the relic abundance 
of this particle is its coupling to the Z boson. We will 
assume a coupling (g/2cw)t- The parameter e fixes the 
annihilation rate of \ pairs into quarks and leptons (with 
a Z in the s channel) and then the contribution of x to 
the energy density of the universe. It is a simple exercise 
to obtain that for M x = 70 GeV we must have c ss 0.25. 
Ref. [25I ] provides an explicit realization of this model in 
a partly SUSY scenario. 

We can now estimate the flux of \ pairs produced by 



cosmic rays. Again, the a xx cross section [26| (see the 



Appendix) depends only on the coupling with the Z bo- 
son (we neglect virtual effects of heavier particles as well 
as the production of heavier fields that may decay into x 
plus standard particles). In Fig.[3]we have included <J XX ■ 
The flux of downgoing dark-matter particles would be 
4.6 x 10 -4 pairs per year, squared kilometer and sterad. 
Their energy distribution is shown in Fig. [5l 



IV. NEW PHYSICS IN A SINGLE EVENT OF 
ULTRAHIGH ENERGY 



We would also like to find the probability for the 
production of gluino or WIMP pairs in single events 
of the highest energy (the ones to be measured at 
AUGER) . In Table U we give two simulations of show- 
ers of high (xp ss 0.3) and low (xp w 0.03) elasticity, 
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FIG. 6: Probability to create X = gg, \X by a primary nu- 
cleon (lower lines) or by the primary or any secondary hadron 
inside the shower (upper lines) as a function of the shower en- 
ergy. 



both with a proton primary of energy 10 10 GeV. For 
M = 200 (300) GeV the probability to produce a gluino 
pair in each shower is ~ 5.4 x 10~ 5 (9.7 x 10~ 6 ) and 
V\ g w 4.3xl0~ 5 (6.6xl0~ 6 ), respectively Around a 20% 
of this probability corresponds to the primary proton, 
whereas the rest expresses the probability that the gluino 
pair is produced by secondary hadrons inside the shower. 
Notice that the shower with no leading hadron after the 
first interaction has more hadrons of intermediate and 
lower energies. For a light gluino these tend to compen- 
sate the effect of the leading hadron (i.e., V~ g ~ T^sg)- 
However, as the gluino mass M increases the less ener- 
getic hadrons are not able to produce gluinos and 
becomes significantly smaller than V~g. 

The probability to produce the pair of WIMPs in these 



two showers is VZ X 



2.5 x 10~ 8 and V b xx » 3.1 x 10" 



In this case the exotic particles, of mass M x = 70 GeV, 
can be effectively produced by hadrons of intermediate 
energy and V xx > V xx . 

In Fig. [6] we plot the average probability to pro- 
duce gluinos (and WIMPs) for a primary proton of en- 
ergy above 10 8 GeV. We separate the probability that 
the gluinos are produced by the primary or secondary 
hadrons inside the shower. We find that secondary pions 
are the main source of gluinos in showers above 10 8 GeV. 
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V. DETECTABILITY OF THE GLUINO PAIRS 

The small flux of WIMPs that we have obtained does 
not seem to imply any experimental consequences in neu- 
trino telescopes or air shower experiments. The flux of 
gluino pairs, however, might be visible. 

As mentioned above, IceCube (around 1 km 2 of area) 
[l3| would be exposed to one gluino pair per year if M w 
160 GeV. Since this mass is already below the experimen- 
tal limits from the Tevatron, such an event would be un- 
expected. For example, if M — 200 (300) GeV the prob- 
ability to have the two gluinos crossing IceCube in one 
year is just around 0.22 (0.015). Let us assume that such 
unexpected event actually happens. The mean energy of 
the gluino pair is around E = 2.1 x 10 6 (5.1 x 10 6 ) GeV, 
and we obtain an average angle between the two gluinos 
9 SS « 8.8 x 10~ 4 (5.4 x 10~ 4 ) rad. We will approximate 
a constant interaction length in water A^, w 175 g/cm 2 
for the gluino hadron and a linear energy loss per inter- 
action AE = kj with k w 0.2 GeV. This means that a 
G of M = 200 GeV deposits a 1 per mille of its energy 
every 1.8 meters of water. It is then easy to deduce that 
the energy of the gluino pair at a depth x is just 



10" 



E. 



G ~ E a exp{-a;/a:o} , 



(6) 



where £ = A^M/fc » 1.8xl0 5 (2.6xl0 5 ) g/cm 2 for M = 
200 (300) GeV. Therefore, two vertical gluinos would 
reach a depth x\ — 1.4 x 10 5 g/cm 2 (the top of IceCube) 
with a total energy around Si = l.lxlO 6 (2.1xl0 6 ) GeV. 
If they are produced at an altitude H = 20 km, the 
gluinos will enter IceCube separated by 18 (11) meters 
and will deposit AE = 4.7x 10 5 (1.5 x 10 5 ) GeV of energy 
through the kilometer long detector. Such a signal could 
be clearly distinguished from a typical muon bundle in 
an air shower core. 

More frequent events could be obtained in the larger air 
shower experiment AUGER. The number of gluino pairs 
of M = 200 (300) GeV hitting per year the 3000 km 2 
AUGER area would be around 



dcostf 2irA cos 9 T $ 



99 



330 (22) , (7) 



where 9 is the zenith angle. However, this is not a very 
significant number, as these gluinos are inside the air 
shower where they were produced, and many of these 
showers have an energy below the 10 8 GeV threshold in 
AUGER. In Fig. we plot the distribution of the 330 
(22) gluino pairs as a function of the primary energy. We 
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FIG. 7: Distribution of the gluino pairs as a function of the 
shower energy. 



obtain an estimate of 20 (2) gluino events per year in- 
side showers above 10 8 GeV. A 75% of these events come 
from a zenith angle below 60° whereas the remaining 25% 
are showers reaching AUGER quasi-horizontally, with a 
larger zenith angle. 

Let us then analyze a typical gluino event inside a 
10 10 GeV shower for M = 200 GeV. Once created, the av- 
erage angle between the gluinos is 4.7 x 10~ 4 rad, whereas 
the total energy of the pair is around 5.2 x 10 7 GeV. The 
gluino-gluon R-hadrons will interact every 160 g/cm 2 , 
depositing in each interaction an average one per mille 
of their energy (this corresponds to the value k = 0.2 
given above) . Notice that the signal produced by a gluino 
hadron is quite different from that produced by a light 

A proton deposits 



27, 2 



hadron of the same energy 
most of its energy within a couple of vertical atmospheres 
(2000 g/cm 2 ), whereas the gluino would interact in the 
same interval around 12 times, depositing a fraction 10~ 3 
of hadronic energy every time. The gluino signal is more 
homogeneous, after 4 or 5 interaction lengths it is a trace 
of constant energy. Notice also that the number of unions 
from pion decays that the gluino produces grows propor- 
tional to the length of that trace. This will affect the cur- 
vature in the shower front hitting the ground. The series 
of mini-showers started by the gluinos will produce muon 
bundles at different slant-depths, and the front will ex- 
hibit a more pronounced curvature than what one would 
expect for a shower of the same energy interacting only 
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in the upper atmosphere. For isolated gluino showers 
(i.e., not inside a proton event) these features have been 
discussed by Anchordoqui et al. 27J and confirmed with 
a Montecarlo simulation in [2^ . Previous works on air 
showers started by (lighter) long-lived gluinos were pro- 
posed as a possible explanation of the events beyond the 
GZK cutoff 0. 

Therefore, one needs to distinguish two types of events, 
the vertical ones (with a zenith angle 8 < 60°) and the 
quasi-horizontal ones (with 9 > 60°). If the gluinos are 
produced in a 10 10 GeV vertical shower they will take a 
0.5% of the initial energy and deposit a 6% of that en- 
ergy in two traces separated by « 9 m. Each trace is the 
superposition of hadronic showers of 2.6 x 10 4 GeV sepa- 
rated by 160 g/cm 2 , i.e., around 1.6 x 10 5 GeV of energy 
in each gluino trace. It seems unlikely that such a distor- 
tion in the profile of a vertical shower can be detected in 
AUGER. Quasi-horizontal showers, however, could pro- 
vide a more promising signal. At large zenith angles the 
depth of the atmosphere increases up to 36000 g/cm 2 for 
6 = it/2, and this has several important effects on the 
gluino signal. First, the two gluino hadrons will cross 
the AUGER surface with a larger separation than in ver- 
tical showers. In particular, horizontal events come from 
a distance D w \/2HRt, where H « 20 km and Rt is 
the radius of the Earth. This gives a separation w Ddgg- 
Second, after a certain depth most of the hadronic energy 
of the shower has been absorbed by the atmosphere, leav- 
ing only muons (and invisible neutrinos). In contrast, the 
hadron content in the gluino trace is basically constant, 
as each G starts a 2.6 x 10 4 GeV shower every 160 g/cm 2 . 
Therefore, the detectors would observe charged hadrons 
produced by the two gluinos in addition to the muons. 
The third effect on these inclined events is that the total 
number of muons produced in the gluino traces is larger 
than in vertical showers, as it grows almost linearly with 
the depth. In this way, the two gluinos inside the 10 10 
GeV event would deposit around 10 7 GeV (a 20% of their 
total energy) if they come horizontally crossing a depth 
of 36000 g/cm 2 . Notice that these muons produced deep 
in the atmosphere increase the curvature of the shower 
front. 



VI. SUMMARY AND DISCUSSION 



accelerator experiments. The primary nucleons or sec- 
ondary hadrons inside an air shower may produce new 
massive particles when they interact with atmospheric 
nucleons. We first have determined the total flux of 
hadrons within cosmic rays. We think Fig.[T]is a relevant 
result, as it is this (and not just the flux of primaries) 
the relevant flux to evaluate the production rate of new 
physics by cosmic rays. Then we have discussed how 
to calculate the creation of a massive particle, focussing 
on the long-lived gluino of split-SUSY models. We have 
shown that the flux of gluino hadrons produced by very 
energetic cosmic rays allows the possibility of isolated 
gluino events crossing a km 2 neutrino telescope like Ice- 
Cube. In the larger detector AUGER, for a gluino mass 
M = 200 GeV we obtain an estimate of 20 gluino events 
per year inside air showers of energy above 10 s GeV. 
We have argued that the gluinos could modify substan- 
tially the profile of the showers with a large zenith an- 
gle. The gluino hadrons propagate losing a 1 per mille 
of their energy in every interaction length, defining an 
air shower with an approximately constant number of 
hadrons and with a number of muons that grows linearly 
with the depth. Therefore, the signal is stronger if they 
are created inside showers that reach AUGER from a 
large zenith angle. We think that the possibility of de- 
tecting in AUGER such a signal deserves a more detailed 
study. 
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APPENDIX A: PARTONIC CROSS SECTIONS 

Making use of the Mandelstam variables, s + t + u = 
2M 2 , and the polar angle 9 in the center of mass frame, 
t = -(1 + P 2 - 2/3cos(9)s/4, with $ = y/l -AM 2 /s and 
7 = (1 — /3 2 ) -1 / 2 , the cross-sections read: 



Cosmic ray experiments could provide a window to ex- 
plore physics beyond the standard model, complementing 



da s(3 da 
dcosfl ~ ~2~"d7' 



° = 2 



dcos# 



da 
dcos6* 



(Al) 
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99 -> 99- gg ->■ XX- 

da 9na 2 s ( 2(t — M 2 )(u — M 2 ) M 2 (s-AM 2 ) 



dt As 2 \ s 2 {t-M 2 ){u-M 2 ) 

(t - M 2 ){u - M 2 ) - 2M 2 (t + M 2 ) 
+ (t - M 2 ) 2 

(t - M 2 )(u - M 2 ) + M 2 (u - t) 
+ s{t - M 2 ) 

{u-M 2 ){t-M 2 )-2M 2 (u + M 2 ) da na 2 e 2 sf3 3 



(u - M 2 ) 2 d cos 6 32 S ^c% (s - M 2 ) 2 



(A6) 



(u-M 2 )(t-M 2 ) + M 2 (t-u) \ x(l-4r ? + 87 7 2 )(l + cos 2 0), (A7) 



S ( u - M2 ) ) '"' ira 2 e 2 sf3 3 M A o 2 , . , . 



4 S lH 1 + 7 _ ^J l0g l-^ 



4 + ^) ; r (A; - !) 



d(j 8^ 2 (t - M 2 ) 2 + (u - M 2 ) 2 + 2M 2 s /4 . 



. with r) = \Qq\sly and e = 0.25 (see discussion in 

— + (A5) Sect. III). 

9s V 7/ 
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